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Abstract 


YVe  analyze  the  long  time  behaviour  of  fully  discrete  solutions  to  a  one-dimensional 
nonlinear  viscoelastic  problem.  It  is  shown  that  these  approximations  which  are 
found  by  a  continuous  time  Galerkin  method  converge  to  a  steady  state.  The  possi¬ 
ble  numerical  steady  states  are  characterized  and  in  particular  their  high  degree  of 
dependence  on  initial  data  and  mesh  design  is  explained.  Computational  results  are 
included  which  show  the  above  dependence  and  indicate  that  the  numerical  solutions 
will  typically  not  converge  to  unstable  states.  -  ! 
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1.  Introduction 


Recently,  attempts  have  been  made  at  minimizing  nonconvex  stored  energy  function¬ 
als  by  means  of  studying  the  steady  states  of  auxiliary  time  dependent  problems.  In 
parts  of  the  literature,  the  class  of  methods  is  termed  Dynamic  Relaxation.  This 
technique  involves  solving  a  time  dependent  P.D.E.  and  marching  to  a  steady  state 
which  hopefully  is  a  near  minimizes  'A’e  refer  to  (lj,  (2],  [3 j ,  j 4 ] ,  [5],  Jfi)  and  the 
references  there.  In  this  paper  we  will  study  the  long  term  behaviour  of  certain  fully 
discrete  solutions  to 

Vtt  -  [o(Ux)  +  in  (0,  1)  x  (0, oo) 

V{  o,0-o,  t>0 

(o(Ux)  +  Uzl)(\,t)  or  C  M)=0,  I>0 
U{z,  0)  =  Uq(x).  Ut(x,  0  ~  Vq(.t)  in  (0,1) 

which  commands  some  interest  in  its  own  right.  (1.1)  models  the  one-dimensional 
motion  under  zero  body  forces  of  a  nonlinear  viscoelastic  material  of  rate  type,  some¬ 
times  called  the  Kelvin- Voigt  model.  U(i.t)  denotes  the  displacement  at  time  t  of  a 
particle  having  position  x  in  some  reference  configuration.  See  also  [7j.  It  is  hoped 
then  that  for  t  sufficiently  large,  U  will  approach  a  local  minimum  and  satisfy  the 
Euler-Lagrange  equations  associated  with  the  minimization  problem. 

The  asymptotic  behaviour  of  (1 .1)  was  die  focus  of  two  studies,  (8)  and  (9),  in  which 
it  was  shown  that  weak  solutions  exist  g.obally  provided  Uo  €  Vq  (r.  1/(0, 1) 

and  the  sign  of  the  stresses  for  large  s  is  restricted  as  follows:  a  :  R  — *  R , 

3M  >0,  |  s  \>  M  =>  scr(s )  >  0.  (1.2) 


Furthermore  the  solutions  converge  strongly  to  (local)  equilibrium,  as  t  — ►  oo: 


Ut  — ->  0,  in  Hl(0,  1),  o(Uz)  +  Uzt  — >  0  in  H2( 0,  l) 
Ux(x,t)  —>  Soo(x)  bounded];.'  a.e.,  <7(5^)  =  0  a.e., 


(1.3) 


and  when  a  dynamic  stability  criterion 

~(s)  >  Oq  >  0  for  a(s)  =  0  (1.4) 

as 

is  satisfied,  the  asymptotic  states  are  metastable  in  the  sense  that  if  the  strains  are 
perturbed  by  small  amounts  except  on  a  set  of  measure  e,  then  the  solution  will 
approach  an  equilibrium  with  stiain  equal  ‘.o  the  unperturbed  strain  everywhere  but 
possibly  the  same  exceptional  set.  !/(•,<)  converges  in  the  sense  of  generalized  curves, 
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of.  [  10]  and  [  1 1  ] ,  i.e.  there  exists  a  parametrized  family  of  probability  measures 
{uT}T(r (o.i)  on  R  such  that  g(uz(-,t))  -~l<  isx,g  >  in  L°°(U,l)  for  each  continuous 
g.  Here,  <  u,g(s)  >  denotes  the  value  of  a  probability  measure  u  for  a  continuous 
function  g  on  the  state  domain. 

It  has  been  observed,  e.g.  by  Silling  |3],  that  dynamic  relaxation  as  a  computa¬ 
tional  method  for  a  similar  two-dimensional  problem  may  yield  numerical  solutions 
exhibiting  different  phases  that  can  be  identified  with  appropriate  ranges  of  values  of 
the  deformation  gradient.  In  mixing  regions  between  two  elliptic  phases,  the  gradients 
become  discontinuous  ami  the  solution  was  observed  to  be  highly  mesh  dependent. 

We  are  here  interested  in  studying  (l.l)  when  a  is  not  a  monotone  increasing 
function  so  that  the  stored- energy  functional 

Mv)  =  /  'l>W)dx  (1.5) 

Jo 

is  not  convex.  In  general  there  are  thus  infinitely  many  solutions  to  the  corresponding 
equilibrium  problem: 

{o(u'(x))}'  —  0  (1.6) 

subject  to  mixed,  Dirichlet  and  traction  boundary  conditions 

u(0)  =  /o  and  n(v')(l)  =  P  ( 1 .7a) 

or  Dirichlet  boundary  conditions 

u(0)  =  /o  and  v(l)  =  f\  (1.7b) 

if}  is  a  real-valued  function  in  one  variable:  ifi(v')(x)  =  o(s)ds  and  is  defined 
as  follows.  Let  si  <  s  <  s^,  0  <  a  <  (su  -  si)/ 2  and  Aj,  A2,  A3  >  0.  We  take 
s  —  fi  —  fo  —  P  —  0.  Let  t l>  be  a  double  well: 

ifj(s)  >  0(s/y)  =  ifj{sv)  =  0  nr  *^sL,sUy 
0(-s)  >  A, {s  -  sL)2  fo.  |  5  -  sL  |<  a, 

ip(s)  >  \[(s  -  su)2  for  j  s  -  sv  |<  a,  (1.8) 

0(5)  >  A^2  for  Si  +  a  <  s  <  sy  -  a, 

A3 (s  -  s)2  >  0(s)  >  A2(-s  -  s)2  for  s  £  [sl  -  a,sv  +  a]. 

1$  is  well  defined  on  //’.  Collins,  Kinderlehrer,  and  Luskin  ( 1 2]  noted  that  the  vari¬ 
ational  principle  infue//i  I$(v)  where  u  satisfies  (1.7b)  may  have  nonunique  limits  for 
minimizing  sequences  as  well  as  nonunique  associated  Young  measure  (the  probabil¬ 
ity  measures  mentioned  above),  jll].  The  simplest  limit  deformation,  however,  is  the 
linear  function  /  satisfying  (1.7b).  In  [12]  they  consider  instead  minimizing 

•-4M  “  l+{v)  h  /  (v(x)  -  f{x))2dx 
0 


(1.9) 


Figure  1:  Energy  and  Stress  as  functions  of  Strain 


subject  to  u(0)  =  u(l)  —  0  with  /  the  unique  limit,  deformation  and  a  unique  associ¬ 
ated  Young  measure 

i'r  =  +  (!  -  (1-10) 

where  7  =  (sy  —  s)/(sy  -  sL)  and  6Z  denotes  the  Dirac  delta  distribution  with  sup¬ 
port  at  2.  Note  that  inf//i9u:(i,7)  I^,{v)  =  \n{v6Hi  J^(v).  The  stored  energy  and  the 
associated  stress-strain  law  are  be  depicted  in  fig.  1  (s  =  0) 

Note  that,  although  infue//i  J$(v)  —  0,  there  is  no  minimizer  in  H1 .  In  [12]  it  was 
shown  that  minimizing  J 0  over  a  sequence  of  (uniform  mesli)  finite  element  spaces, 
leads  to  a  minimizing  sequence  «h,  such  that  u!h  converges  in  a  weak  sense  to  the 
unique  Young  measure  in  (1.10)  as  h  -*  0.  u’h(x )  oscillates  (in  the  limit)  between  the 
energy  wells  at  s  —  Sl  and  at  s  =  5/,  occupying  these  states  in  proportions  7  and  1  —7 
of  the  interval,  respectively.  Double  wells  of  unequal  heights  can  be  incorporated  as 
in  (12)  by  shifting  0  by  some  linear  function. 

Example  1.1  Let  us  characterize,  for  future  use,  the  solution  to  the  discrete  Galerkin, 
Euler-Lagrange  equations  in  the  simplest  case.  Let  m  £  N,  for  h  =  ,  i0  —  0  and 

x,  =  i/i,  /,  =  (xj_ j,xt)  for  t  —  l,  -,m.  Define  the  finite  element  space  for  the 
boundary  conditions  (l.7a)  (f0  =  s  —  0,  I  —  (0, 1)) 

S]  =  C(0,1)  :  u(0)  =  0  and  v  |/.(E  Pi (/,•),  t  =  l,---,m}  ( 1 . 1 1 ) 

Then  the  weak  form  of  the  Euler-Lagrange  equation  (o^u*)),,  =  0  becomes 

Find  Uh  £  Slh  such  that  for  all  /  a(uh  z)vx  dx  =  0.  (1.12) 

Since  i/h(0)  =  0  and  v.hiZ,  vz  are  piecewise  constants,  we  can  substitute  £),•  o(uh,z)(v{x,)~ 
i'(x,.  ,))  for  /,.  Now  successively  testing  a.gainst  rn  basis  functions  ?V  €  Sjx  defined  by 
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v'(xt)  --  1,  v'(tj)  ~  0  for  j  i,  i  —  1  ,  one  easily  see  that  o{u'h )  remains  the 

same  constant  throughout  I .  So  u/,  £  C]h  where 

C'h  {»'  f  :  o(vz)  ~  0}  -  {t-  €  Slh  :  vz  £  { sL,0,sv }},  (1.13) 

such  that  Cl  consists  of  3'’1  isolated  solutions.  Any  two  nonidentical  solutions  in  C\ 
differ  by  at  least  a  positive  constant  depending  on  h,  s i,  and  sp. 


Example  1.2  In  the  case  of  Oi  rich  let  boundary  conditions  at  both  endpoints  (1.7b), 
the  discrete  equilibria  become  less  tangible.  With  m,  h,  x,  and  /,  as  in  the  previous 
example,  define 


sh,0  ~  (r  e  <?(°>  *)  :  l'(°)  ~  y(0  =  0,  v  |/,.e  Pi  (/,•),  I  =  i,  •  •  ■  ,  m}  (l.  14) 

as  above  we  see  that  a(uhz)  is  constant  P,  say.  Let  o~x[P)  =  {yU],/r2,/i  3}  listed  in 
ascending  order.  Let  m,  denote  the  numbers  of  intervals  in  which  u^  T  —  p,-,  1  —  1,2,3. 
Since  u ^  satisfies  (1.7b)  we  get  the  following  constraint 

3 

]>_)  m,g,  =  0,  m,  >  0  for  i=  1,2,3. 

t=i 

The  set  of  possible  values  for  P  depends  very  much  on  the  shape  of  the  stress  strain 
law.  It  is  clear  that  o~l(P)  =  a  singleton  is  not  viable,  similarly  if  a_1(P)  = 
(p-i,  P2},  we  need  pi/p2  £  Q  and  m  sufficiently  large.  To  explain  take  a  very  simple 
minded  trilinear  stress-strain  law,  such  that  a  has  roots  sp,  0,  and  sn\  j  o'  |=  A.  Thus 


c(s) 


A(s  -  sh),  for  s  <  sij 2 
-As,  for  si/2  <  s  <  sp/ 2 

A(s  -  sp),  for  s  >  sp/ 2 


In  one  of  the  two  root  cases,  P  =  — Asp/2  and  n  =  sp/2,  p2  —  Sp  -  sp/2  so  that 
Si  j  su  must  be  rational  (-^,  irreducibly  say)  and  m  >  p  +  q.  In  the  three  root 
cases  -P/A  £  |(-.sp,-sL)  and  p.)  =  sL  +  P/A,  p2  —  -P/A,  and  p3  =  sp  +  P/A. 
the  constraint  of  zero  mean  slope  above  can  be  met  always  by  letting  m2  =  0  and 
choosing  P/A  =  -(m(Sf,  +  m3.su)/m.  This  allows  for  at  least  one  degree  of  freedom 
in  choice  for  P  and  thus  possible  slopes  p,,  as  long  as  m,  >  m/4,  t  =  1,3.  It  is 
particularly  interesting  that  we  can  allow  m2  >  0  putting  strict  lower  bounds  on 
m3.  The  solutions  are  still  isolated,  the  P/A  values  admissible  being  separated  by 
dh  =  c($i,su)/m  amounting  to  a  nearest  neighbour  distance  of  the  same  order  in  the 
set  of  discrete  Galerkin  equilibria.  Of  course,  more  nontrivial  discrete  equilibria  are 
realizable  using  nonuniform  meshes  selected  appropriately. 

Increasing  the  polynomial  degree  in  the  finite  element  space  will  likely  yield  an 
even  wider  spectrum  of  discrete  equilibria. 


5 


Other  forms  of  damping/dissipation  could  also  be  considered  for  dynamic  relax¬ 
ation:  thermal  (heat  diffusion),  frictional  (including  Ut  in  ( 1 . 1 ) ) ,  and  viscoelastic 
of  history  type  ( o{Ux )  depending  on  the  deformation  for  all  previous  times  through 
convolution  with  a  kernel  (see  Bielak  and  MarCamy  [13]).  Dissipation  mechanisms 
which  will  substantially  influence  behaviour  include  capillarity  (adding  a  Uzxxx  term 
to  (1.1),  see  Slemrod  [  15])  and  non-local  (in  space)  constitutive  relations  ,  see  Be- 
lytschko  and  Bazant  [5).  These  are  typically  not  used  in  dynamic  relaxation  with  the 
exception  of  Beiytschko. 

Finally  note,  that  if  the  energy  in  contrast  is  convex,  a  unique  solution  is  given  by 
the  Kuler-Langrange  equations  and  there  exist  numerical  procedures  for  minimization, 
see  e.g.  [16],  and  the  time  dependent  P.D.F.  was  handled,  see  [17],  [  1 8] ,  [19].  If  one 
uses  i, he  Maxwell  relation  to  relax  the  problems,  nonuniqueness  of  an  even  larger  class 
than  before  arises.  Numerical  methods  exist,  cf.  (20 j  and  the  references  there. 

Our  goal  will  be  to  establish  the  set  of  conditions  under  which  a  class  of  numerical 
methods  will  yield  a  long  time  behaviour  with  asymptotic  states  that  are  (local) 
minima  of  (1.9). 

The  plan  of  the  paper  is  as  follows.  We  analyzed  what  possible  numerical  steady 
states  exist  for  (l.l)  already.  In  the  next  section,  we  introduce  the  continuous  time 
Galerkin  scheme.  In  section  3,  we  show  that  the  numerical  solution  must  converge 
to  one  of  the  steady  states  (for  the  specific  fully  discrete  method  in  section  2).  We 
report  on  some  of  our  numerical  experiments  in  section  4  which  is  followed  by  some 
concluding  remarks. 
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2.  Dynamic  Relaxation  by  Continuous  Time  Galerkin  Schemes 
with  Viscoelasticity 

We  introduce  the  continuous  time  Galerkin  (CTG)  schemes  (see  [22]  and  the  refer¬ 
ences  there).  Let  0  =  x0  <  .r,  <  ...  <  x,  ,  <  x,  <  ...  <  xm  =  1,  /,  =  (x,..1(  x,). 
hi  1  /,  |,  i  ---  1 ,  ■  •  ■  ,  m  and 

s'h  ---  (x  G  c(0,  1)  :  x(0)  -  0,  x  |/,-€  P\{P),  i=  b  •  ■  •  ,m}. 

Let  0  —  to  G  1 2  • .  • ,  •/„  (tn—  i >  Li ) i  bi  |  j i  ^  0,1,2,...  and 

•h'fc  =  {t  G  C(0,oo)  :  r  |j„€  /’i(d.,)}. 

Then  5^  -  Slh®S'k  -  (A  :  A  --  )T'  =  |  X;b>  X;  G  5/(,  G  .9*,  l  6  /V}  and  we  discretize 
as  follows  (CTG): 

Find  u  6  S/lJt  and  v  6  .9^  such  that 

(K  -  *>,  x))n  =  0  Vx  €  Si  ®  Po(^n) 

(2.!) 

((uo  A))„  +  ((u*j  +  <7(ux),  Ax))n  =  0  VA  E  5”^  ®  Po{Jn) 

where  ux(-,  0)  -  u0>1  =  U0  x  and  u(-,  0)  =  v0  =  V'n. 

The  inner  products  are  defined  by 

(v,  w)  =  fj  u(x)w(x)  dx 


[v,w)n  =  fJn  v(t)w(t)  dt  (2.2) 

((u,  w))n  =  fo  fjn  v(x>  t)w(x,  t )  dtdx. 


Note  that  (2.1)  can  be  given  the  following  equivalent  finite  difference  in  time  formu¬ 
lation 

/o(u"+i"“"  -  ^'r")xdx^0i  Vx  G  SI 


0  \  k 


x  +  dx,  VA  e  si. 


where  w}  —  w(tj)  for  j  —  n,  n  4 -  1;  w  =  u,  ux,  v  6  S^.  One  of  the  fundamental 
properties  is  the  following  energy  estimate. 
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Proposition  2.1  Let  u,  v  t  Slhk  be  the  solutions  of  (2.1).  Then  for  any  n  6  N,  the 
following  energy  identity  holds 

f  (,!”2  +  d>{uz))dx  -  /  /  ul'dtdx  (2.1) 

■'o  2  Jo  • 


Proof :  Choose  \  •-  vt  and  A  -  ut  in  (2.1)  to  got 

((«t,  v())n  -  l  J  V2dx  i;-;-', 

2  J  o 

and 

((tp,  v<))„  =  -((^K)  »  »io  «i())n  =  ~  xl>{ut)dx  | 

•'o 

from  which  ('2.1)  follows. 


<» 


o  ^ 


iij,  dtdx. 


□ 


Corollary  2.2  Under  the  same  hypothesis  as  in  Prop.  2.1, 

X  fi 


E(tn)  -  Eo  -  /  I  u]tdxdt 
'o  •'o 


(2.5) 


where 

E{t,i)  -  fo{\v2  1-  0(7xI))(.r,  ftl)  r/x, 
Eo  E(t0). 

Proof :  Sum  (2.1)  from  j =  0  to  j  =  n  —  1. 


(2.6) 


□ 


We  still  have  to  prove  existence  and  uniqueness  of  solutions  to  (2.1) 

Lemma  2.3  If  (u,u)  if  a  solution  to  (2.1),  then  there  exists  an  M  >  0,  depending 
only  on  Eq,  such  that  for  all  n  C  N 


(•)  IKIL’io.i)  <  m 

(•*)  lltWOlU-,,.,,  <  m 

(lit)  ([tt,‘||//i(o,t)  <  M 


(2.7) 
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Proof-.  (2.5)  yields  jUi"||^(o.i)  <  (2 E1,  2)  and  ||d>(u")IU'(o,i)  -  from  winch  («) 
and  (n)  follow.  Using*  (it)  and  subdividing  /  according  to  which  of  the  last  four 
sets  defined  in  (1.8).  m"  belongs,  the  I.2  norm  of  u"  is  bounded  and  (in)  follows  by 
Poincare’s  inequality. 

1  : 

We  then  phrase  (2.1)  as  a  fixed  point  problem.  Consider  solving  the  CT(!  problem 
on  the  time  slab  Qn  -  (0,  1)  x  ./„.  Define  the  map  :  (Sf  <g>  P\  ( J„))2  ( S]h  2 

as  the  solution  (u; * 1 .  v1 1 1 ),  given  (tP.r’),  to 

Find  (-  ,V‘  X)  Pi(J„)  such  that 

((ur1  -  \i).  o.  v> .  sih 

(2.8) 

((■r'n))„  -  A,))„  -•  -((<>(«'),  A,)).,  VA  6  S> 

wliere  u2+I  [x,  t)  =  u"(x)  and  iP  +  1  (x,  f„)  =  vn(x),  Vx6(0,l). 

This  defines  a  linear  system  of  equations.  To  show  there  is  a  unique  solution  on  each 
iteration  we  prove  the  homogeneous  system  (zero  right-hand-side  and  zero  data  at 
t  =  tn)  has  only  the  zero  solution. 

Lemma  2.4  Let  <l>  be  defined  as  m  (2.?).  Then  <I>  1  (0)  -=  0. 

Proof :  With  zero  right-hand-side  in  (2.8),  let  A  —  u}t  *  1  and  x  •-  v{ 1 1  so  that 

-  /  (v;  +  1(x,tri  t  l)):d£  !  /  [  u3zl  '  dr.dt  -0 

2  Jr>  J  J„  Jo 

and  iiJxl 1  =  0,  u;  +  1  =:•  0,  and  fd+ 1  (<n+i )  =  0.  Now  choose  A  —  vj  hl  to  get  ((f;/41 ,  v{ ' 1 )),, 
0  and  v{+l  =  0  and  v,n  =  0. 


Bn  -  {x  €  Si®  Pi(Jn)  ■  max||xx(-,OII/>2(o.i)  ~ 

and  Mc  =  sup|,|<M  +  1  |  ct(s)  |.  Then  $  selfmaps  B ^  into  for  R  —  M  +  1,  M  being 
the  constant  in  (2.7), 

Lemma  2.5  let  <I>  be  defined  as  in  (2.8).  Then  <f>  :  B’f^  ,  i  -  flJEr.,,  provided  kn( A-/2  -I- 
knMf)  <  1. 
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I'roof :  As  in  I, ho  previous  proof,  let  A  tij f  1  and  y  -  'V  * '  to  get 

'  f  (vJ")J(/:r  li;;*'  •  (('0/ '))„  -((<r(uJr),uithl))n 

L  ~  it 


•  1  \2I  t  .  \,/v  .  /  1  I 


(iV  )‘  |  j.  f .  i  )r/:r  -  /  J  (iVT,  )  dt dr 

(i 

/  (i‘)‘dr  l  /  j  (o(x~r [))2dxdt 

•'ll  •'it  ■'j,, 

•  A/2  •  a-„a/; 


l|w;x"(-.0  «','!/.•  1  ■  *r(/  J  (K'fdtdx)'1' 

•II 


-i)  ’  hat 


ll«i’  l(-.OII/.a(r..D  •  K  /.  .1  *  {A-n(A'/2  h  knXll))xl2  <  M  |  1 

ending  the  proof. 


Finally,  <P  is  a  contraction  in  some  appropriate  norm.  Let  La  -  sup|,|<AY.( ,  |  or' ( .s ) 
and  l!|u|||2-  {(vr,vT))„. 

Lemma  2.0  Let.  <f>  be  defined  as  r  fC.  *).  Then  I  k  <_  (0,  l)  such  that 


l!l<M 

f/|.  ?-,) 

<!’ :  1  n  2-  ''2)111 

'  K  |f|(M| 

U2||| 

(2.9 

and 

|  [  (F  2  (  u  1  ,  t  ’  I 

J  ^2  (" : 

•  )  (0,  •  i)|| 

/-2<0,1)  k 

lll«i  -•  u2|l| 

provided  kj/2 [jt 

n  <  -  1 1  'tf  ( n , , 

<’,)  •  in, 

t  ! ,  2 

I’roof  :  Let  6 

-m  ~jw 

U  ,  U  2 

and  f 

---I  -j  *  1 

'"1  ’h  ■ 

Then  from  (2.8), 

(2.8)2, 

we  get 

((*i  ' 

•  X  )  )u  --=  0 

VX  €  5* 

Uo 

i  A))rt  d  (  (hjt 

,**))„  • 

a(u;2  J, 

A*))„  VA 

6  S* 

Letting  v  ((  and  A  f>(,  we  get 


~  llf  (■’  Lx  1 1 1  /.-’(<'  i)+  I|| Sz(  ||j 

<  Ni,  -  “Ul  infill 

■  \  i1  *  2  ^  "z.xlll2 
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Since  !l|tfrj||<  /c,1/2  |||<5I(||!,  we  get 

l#MI<  llluj,  -  nun2 

yielding  (2.9). 


This  establishes  the  existence  and  uniqueness  of  the  solution  to  (2.J)  by  applying 
the  contraction  mapping  principle  to  (2.8). 
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3.  Large  Time  Behaviour  of  Numerical  Solution 

We  now  analyze  the  behaviour  in  the  large  of  solutions  to  Continuous  Time  Galerkin 
numerical  method  introduced  in  Section  2. 

If  we  again  denote  un  —  u(t„)  (and  similarly  for  vn),  we  next  state  and  prove  the 
following  properties  of  the  sequence  {(u",v")}. 

Proposition  3.1  Let  ( un,vn )  be  the  sequence  uniquely  determined  by  (2.S)  and  u 
the  corresponding  discrete  function.  Then  3 M  >  0,  such  that 

(0  lkr'lk-(o,i)  <  Vn  6  N 
(*0  110 «)ik '(0.1)  <  A/,  Vn  €  TV 

(Hi)  ||tin!k'(o,i)  <  Vn  €  N  (3.1) 

(iv)  /0°°  /0‘  |  uxt  |2  dxdt  <  E0 

(v)  SJn  fg  !  uzl  \2  dxdt  —  *  0  as  n  — >  oo. 

Proof :  (t)  through  (Hi)  merely  restate  lemma  2.3.  Using  (2.5)  in  the  limit  n  — »  oo 
yields  (tv)  which  has  (v)  as  an  easy  consequence. 

□ 
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whore  An  =  {V1  f  vn  1 } /2 .  Note  that  un+1  —  vn  1  =  2(An+l  -  An).  Summing  (2.1)2 
we  have 

H(Wux),  A))fl  +  ((<7(us),  A))rl_! 

(  (  )  )n  H“  (  {^zt  ?  ^ x  )  ) n-  1  | 


+  /'a 

J  r\ 

r>+! 


•.i  +  l 


i„- i 


u(  d£)  dx 


n+l  f  i 

£((«,£,*,)),- +  4  J  (A"n  ~A")*dx. 


;  =  n 


From  (3.2. i)  and  (3.3)  we  have 


n+l 


lim  ]C((a(u*)'  Ax))i  =  °- 


)-n 


But 


uz(t)  =  u"  +  uzt(t  -  tn )  for  t  €  Jn+ 1  and 
uz(t)  =  K  +  uz‘(tn  ~  t)  for  t  6  Jn. 


Thus 

where 

So 


a{ux)  =  a(u")  +  R 

11^11^(0,1)  <  Lak\\uzt\\L7{0}l)  -»  0  as  n  -»  oo. 

n  +  1 

Z^(((7(U*)' Ax));  =  2A:(a(u"),As)  +  i? 

>=« 

where  /?  — »  0  as  n  — ♦  oo.  Thus  (*"*)  holds. 


□ 


Thus,  if  {u"}  or  a  subsequence  tends  to  a  limiting  function,  the  latter  must  belong 
to  Cj[.  To  see  that  (u'1}  converges,  we  proceed  in  two  steps. 

Proposition  3.3  Let  {un}  be  the  sequence  uniquely  determined  by  (2.8).  Then  there 
exists  u  €  C l  and  p  >  0  depending  on  k  the  time  step  such  that 

||un+1  —  u * ||  <  M||u"  —  u" ||  (3-3) 

where  ||u||  =  ||u||/^(0,i)- 
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Proof :  Let  [un,vn)  denote  the  solution  to  (2.3)  and  let 


o : 


O’) 


(3-4) 


for  some  u'  6  Cf  to  be  determined.  From  (2.3) 

0(«;u  l)  -  '!•<.<) 


1C1.*  )  +  q 


u 


n+  i 


-i-  «r'.A,)  =  (»:,A)  +  K,A,) 


so  that 

(C1^)  +  (~kD  +  0Z\\x)  =  {0”,  A)  -  (0UiX,  Xz)  (3.5) 

where  we  denote  D  —  o{u'T)  -  (t/>(u£+I)  -  t/>(ur))/(ux+1  -  u").  Note  that  u  6  Cf 
satisfies  the  discrete  Euler-Lagrange  equation.  Choosing  A  =  0"+1  in  (3.5)  yields 


(orix+1)  -  wx+l)  +  lie  n2  =  (^.c1)  +  (c-c1)- 


Likewise 


( 


u 


n+1  _  ,/n 


I 

u"  v  ,vn+l  +  V 


(3.6) 


n+l  I  n 

»X)  =  (- - r - >x) 


which  with  x  ~  t'n+1  ~  vn  becomes 

0n+l  _  0n 


nn+l  _  Qn  1 

(  -  K)  =  i((C)’  -  W)2.  I) 


or 


(C.o  -  (cc1)  =  jiir'ii’  +  («c.«:,+1)  -  <ca  -  jin 


which  we  substitute  into  (3.6)  and  get 


jiicii’  +  iicir  = 


mac1)  +  O 

(3.7) 

-ice1) +  w.a  + *  iwiii’ 

(3.8) 

Expanding  i/>(u"+1)  about  uf  in  a  Taylor  series;  using  the  previously  established 
bounds  on  u”,u”+l,  as  well  as  noting  |u*|  <  max{|s£,|,  |^t/|}i  and  assuming  o'  is 
Lipschitz  gives  the  following  inequality: 


\D\  <  C‘(IC.xl  +  lCx‘1) 


14 


where  C  is  a  constant.  Applying  the  Schwarz  inequality,  the  elementary  inequality 
(2  xy  <  ex2  +  r  iy2,  x,y,e  >  0),  and  Poincare’s  inequality,  we  can  absorb  terms 
involving  superscripts  n  +  1  on  the  right  hand  side  of  (3.8)  into  similar  terms  on  the 
left  and  we  arrive  at 


Mior + 3ii*ru’  <  ,‘(k)  +  3ii«;i!!),  v>.  e  n 


which  is  (3.3)  and  1  <  /.i(k)  <  Ck  independent  of  n. 


□ 


Theorem  3.4  Let  k  <  k‘  and  (u",  vn)  be  the  uniquely  determined  sequences  from 
(2.S).  1'hen  3u'  £  C l, 

lim  (un,un)  =  (u\0)  (3.9) 

n  -♦oo 

Proof :  We  adapt  an  argument  in  [23].  From  (3.1.m),  || tz" j| // » (0,x)  is  bounded  and  we 
may  extract  a  subsequence  {un>}  converging  to  some  u‘  £  S l-  Because  of  (3.2.i), 
||url+1  —  un||  — >  0  as  n  — ►  oo.  Then  also  u"'  +  1  — *•  u*  as  j  — ►  oo.  Then  (t/>(u"j  +  1)  - 
4>(unP))/(unx’  +  l  -  u"')  — >  0  as  j  — >  oo  and  (er(u).),  A*)  =  0  following  (3.2. u)  so  that 
u*  £  Cl  and  has  to  be  one  of  the  isolated  points  there.  Let 

B{u\e)  -  {x  £  Si  :  ||x  -  u*||  <  (} 

and  pick  S  <  dhf  2  (cf.  Example  1.1)  so  that  u*  is  the  unique  element  in  ClnB{uw  ,26). 
Let  {ur9}  be  a  subsequence  of  {u'1}  such  that  {u">}  C  B{u' ,6 / n).  Then 

{un>+1}  C  B{u,6) 

according  to  (3.3).  We  then  distinguish  between  two  mutually  exclusive  cases: 

Case  1:  3{u"*}jt6/v  C  B(u' ,6)\B(u~  ,6/ n).  Then  this  new  subsequence  must 

have  a  limit  in  B (u* , 6)\B (u* , <5 / yf)  which  contradicts  that  u*  is  the  only  element 
in  B{u\26)nCl. 

Case  2:  ->(3 [unk}keN  C  B[u* ,6)\B{u'  ,6/ fj,)).  Then  the  entire  sequence  { un } 

converges  to  u\ 

□ 
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4.  Numerical  Results 


Jn  this  section  we  compute  the  solutions  obtained  by  the  CTG  method  in  order  to 
clarify  our  results,  demonstrate  the  instability  present  in  equilibrium  values  /c  where 
o' (/i)  <-.  0,  and  show  that  the  large  time  solutions  are  dependent  on  both  the  mesh 
and  the  initial  data. 

To  implement  the  CTG  scheme  we  eliminated  the  velocities  from  (2.3)  to  obtain 
the  following  three-step  method:  Find  C  Si  such  that 


(*«,<+') +*«' 


u"))  +  dtunz]Xx)dx  =0 


VA  6 


S 


i 

h  i 


where 

d\un  =  ^2(u"+1  -  2un  +  u”'1), 

dtun  =  -7(«n+l  -  u"“‘). 

2  ky  ’ 

,  0(u)  -  t p{v)  .  w  , 

o[u,  v)  —  — — - ,  and  rp  (u  j  =  u(u  -  S/Jfu  ~  Sy). 

We  factored  the  numerator  in  6  to  eliminate  u  —  v  from  the  denominator.  We  have 
assumed  uniform  time  steps.  To  solve  this  nonlinear  system  on  each  time  step  we 
used  the  following  fixed  point  iteration  where  we  searched  for  — >  un+1: 

/'(2(y+1)  A  +  ~z[i+^\x)dx 

V  2  1  x> 

f  i  t 

-  /  {(2un  -  un_1)A  +  u"_1A,}rfi 
■'o  2 

+  v /Ws.T)  <))*.«<*  VAe  S> 

2  Jo 

On  each  iteration  the  linear  tridiagonal  system  of  equations  was  solved  by  Gaussian 
elimination.  Since  our  primary  interest  was  in  the  behaviour  of  solutions  of  (2.3)  we 
continued  the  iteration  until  ||z(J  +  1)  —  z^\\lo°  <  10-8  rather  than  carry  out  a  more 
efficient  incomplete  iteration  procedure. 

Figures  2,  3,  and  4  show  large  time  solutions  (which  we  believe  to  be  the  final 
states  for  these  problems).  In  all  these  cases  we  took  the  final  time  as  T  =  40,  the 
number  of  time  steps  as  N  =  200,  and  initial  data,  specified  completely  at  the  nodes, 
as 

u°(xj)  =  ul(xj )  =  10-10  sin(20 xy)  for  j  =  0, 1,  •  •  • ,  m. 

We  took  m  =  9, 10,20  in  Figures  2,  3,  and  4,  respectively.  In  all  three  cases  we  had 


KO 


5  x  10 


10 
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We  conclude  that  each  of  the  large  time  solutions  is  significantly  different  and  none 
are  c lore  to  the  zero  solution  which  does  solve  the  Euler-Lagrange  equation,  hut,  is 
apparently  dynamically  unstable. 

Figures  5,  6,  and  7  show  the  approximation  of  the  large  time  solutions  of  the 
problem  with  Dirichlet  conditions  at  both  endpoints.  In  each  of  these  cases  we  took 
T  ■  80.  .V  —  *100,  and  the  initial  data 

t i°{x)  —  10  ~10  sin(207rx;)  for  j  =  0,  1, . . . ,  m. 

We  took  m  —  9,10,20  in  figures  4,  5,  6,  respectively.  In  addition,  we  found  P  = 
.2078265621  ±  5  x  10~10,  .323969547  ±  5  x  lO”10,  0.0  -  5  x  10'10  in  the  three  cases. 

We  tried  a  larger  number  of  timesteps  in  some  of  the  cases  above  and  found  no 
change  in  the  steady  state  patterns  formed. 

A  possible  explanation  of  the  instability  of  the  zero  solution  was  observed  in  [ 8 J 
by  studying  a  linearized  problem  around  zero.  We  do  the  same,  as  follows:  Find 
w  =  w(x.t)  so 

wlt  =  o'(0)wtt  +  wzzt 

with 

w(0,t)  —  wz(l,t)  =0,  t  >  0 

and 

w(z,0)  -  U0(x),  wt(x,0)  =  F0(x)  in  (0,1). 

The  solution  has  the  form 

OO 

w(x,t)  ~  Y2{AneX"‘  +  Bnex"  ‘)  sin((2rr  —  l)7ri) 

n—  1 

where 

A*  =  1(-(2 n  -  1  )V  ±  yj{2 n  -  1 -  4a'(0)(2n  -  1  )2tt2  ). 

2 

Thus 

eA  +  i  ^  e-c'(0)l 

eK,t  «  eo'(0)-n^  0  as  n  _  0> 

Thus  the  solution  of  the  linearized  system  has  a  component  whicli  is  growing  expo¬ 
nentially. 

A  similar  procedure,  although  more  complicated,  can  be  carried  out  on  the  fully 
discrete  approximation  problem  to  show  the  same  results  when  k  =  ch  where  c  is  a 
constant  and  h  is  sufficiently  small. 
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5.  Concluding  Remarks 

We  have  shown  that  numerical  solutions  to  the  model  problem  (1.7)  found  by  the 
continuous  time  Galerkin  scheme  converge,  as  t  —*  oo,  to  a  numerical  steady  state, 
as  one  would  wish.  If  this  time-dependent  process  is  viewed  as  dynamic  relaxation, 
it  is  succesful  at  attaining  steady  states.  We  prove  and  observe  though,  as  Silling 
[d|  observed  for  a  two  dimensional  anti-plane  shear  problem,  that  these  are  highly 
dependent  on  the  mesh  and  initial  data.  We  can  also  interpret  our  results  as  an 
analysis  of  the  large  time  behaviour  of  the  Kelvin- Voigt  model  for  viscous  damping 
using  a  general,  systematic  approach. 
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